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Abstract 

The Bruggeman formalism provides an estimate of the effective permittivity of 
a particulate composite medium comprising two component mediums. The Brugge- 
man estimate is required to lie within the Wiener bounds and the Hashin-Shtrikman 
bounds. Considering the homogenization of weakly dissipative component mediums 
characterized by relative permittivities with real parts of opposite signs, we show that 
the Bruggeman estimate may not be not physically reasonable when the component 
mediums are weakly dissipative; furthermore, both the Wiener bounds and the Hashin- 
Shtrikman bounds exhibit strong resonances. 
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1 Introduction 

Metamaterials in the form of particulate composite mediums are currently of considerable 
scientific and technological interest [1]. Provided that wavelengths are sufficiently long com- 
pared with the length scales of inhomogeneities, such a metamaterial may be envisaged as a 
homogenized composite medium (HCM), arising from two homogeneous component mediums 
[2, 3]. HCMs with especially interesting properties may be conceptualized if the real parts 
of the relative permittivities (and/or relative permeabilities) of the two component mediums 
have opposite signs [4]. This possibility arises for metal-in-insulator dielectric composites 
[5, 6] and has recently become feasible with the fabrication of dielectric-magnetic materials 
displaying a negative index of refraction in the microwave frequency range [7, 8]. 
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Over many years, several theoretical formalisms have been developed in order to estimate 
the effective constitutive parameters of particulate composite mediums [2]. In particular, the 
Maxwell Garnett and the Bruggeman homogenization formalisms have been widely used [9]. 
Generally, the Maxwell Garnett formalism is seen to hold only for dilute composite mediums 
[10].^ More widely applicable is the Bruggeman formalism that was initially founded on 
the intuition that the total polarization field is zero throughout the HCM [16] . A rigorous 
basis for the Bruggeman formalism is also available, within the framework of the strong- 
permittivity-fluctuation theory (SPFT) [17, 18]. 

Estimates of HCM constitutive parameters generated by homogenization formalisms may 
be required to lie within certain bounds. In particular, the Wiener bounds [19, 20] and the 
Hashin-Shtrikman bounds [21] are often invoked. The Hashin-Shtrikman bounds coincide 
with the constitutive parameter estimates of the Maxwell Garnett homogenization formalism 
[20]. The applicability of theoretical bounds on the HCM permittivity has recently been the 
focus of attention for composites specified by relative permittivities with positive-valued real 
parts [22]. 

In this communication, we consider the application of the Bruggeman formalism, together 
with the Wiener and Hashin-Shtrikman bounds, to isotropic dielectric HCMs which arise 
from component mediums characterized by complex-valued relative permittivities whose real 
parts have opposite signs. This is scenario is typical of metal-in-insulator HCMs [20, 23], for 
example. By duality, our analysis extends to isotropic magnetic HCMs. It also extends to 
isotropic dielectric-magnetic HCMs, because the permeability and the permittivity are then 
independent of each other in the Bruggeman formalism [24] (as also in the Maxwell Garnett 
formalism [25]). Therefore, our findings are very relevant to the application of homogeniza- 
tion formalisms [4] to mediums displaying negative index of refraction [26], for example. 
Furthermore, the implications of our mathematical study extend beyond the Bruggeman 
formalism to the SPFT as well [3]. 

A note on notation: An exp(— 2cjt) time-dependence is implicit in the following sections; 
and the real and imaginary parts of complex-valued quantities are denoted by Re ( • ) and 
Im ( • ), respectively. 

•^The restriction on the applicability of the Maxwell Garnett formalism to dilute composite mediums 
generally emerges from comparison with experimental data [9]. As the particulate volume fraction increases, 
the distribution of particles begins to lose the randomness which is inherent to the theory of the Maxwell 
Garnett formalism [11]. However, the restriction could be bypassed if the distribution of particles in a 
composite medium continues to lack order even under densification, which thought underlies the random 
unit cell approach developed by Smith and colleagues [12, 13, 14]. An anonymous reviewer has suggested 
that self-assembly techniques [15] can yield randomness even at large particulate volume fractions, and could 
therefore extend the applicability of the Maxwell Garnett formalism. 
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2 Theory 



2.1 Bruggeman equation 

Consider the homogenization of two isotropic dielectric component mediums labelled a and b. 

Let their relative permittivities be denoted by and e^, respectively. For later convenience, 

we define ^ 

— if ea,eheR, 

e& 

- , (1) 

The Bruggeman estimate of the HCM relative permittivity, namely e|cM) is provided 
implicitly via the relation [9] 



" /g (66 + 26- J + h (6„ + 26- ' ^ ^ 



wherein /„ and /b are the respective volume fractions of component mediums a and 6, and 
the particles of both component mediums are assumed to be spherical. The Bruggeman 
equation (2) emerges naturally within the SPFT framework [3]. A rearrangement of (2) 
gives the quadratic equation 

2 (e^cM)' + Cm [6g (A - 2/„) + e, (/, - 2^)] - e^e, = 0. (3) 

Only those 6Hc^-solutions of (3) are valid under the principle of causality as encapsulated 
by the Kramers-Kronig relations [28] which conform to the restriction Im (6^^m) — 0- 
Let A be the discriminant of the quadratic equation (3); i.e., 

A = [ea ih - 2Q + 6fe Ua - 2h)f + 86„66. (4) 

Since /& = 1 — /„, we may express (4) as 

A = Qfl (6, - e,f - 6fa {el - Se^e, + 2el) + (e, + 26,,)' . (5) 

An insight into the applicability of the Bruggeman formalism may be gained by considering 
the /(j-roots of the equation A = 0; these are as follows: 



_el- 36„66 + 26^ ± 2V2^-6„66 (6a - e^f 

/g|A=o — 72 • (6) 

3 (6a - 66) 

On restricting attention to nondissipative component mediums (i.e., ea,b € M), it is clear that 
fa\ A=o complex-valued if 5 > 0. Consequently, A > which implies that 6^^^ e M. On 
the other hand, /a|A=o real-valued if 5 < 0. Thus, the Bruggeman estimate 6|cm ^r 
5 <0 may be complex-valued with nonzero imaginary part, even though neither component 
medium is dissipative. 
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2.2 Bounds on the HCM relative permittivity 

Various bounds on the HCM relative permittivity have been developed. Two of the most 
widely used are the Wiener bounds [19, 20] 



and the Hashin-Shtrikman bounds [21] 



(7) 



HSb = ea + 



3/bea(£6-go) 
e(,+2ea-/6(ei,-ea) 



(8) 



While both the Wiener bounds and the Hashin-Shtrikman bounds were originally derived 
for real-valued constitutive parameters, generalizations to complex-valued constitutive pa- 
rameters have been established [23]. 

The Hashin-Shtrikman bound HSq is equivalent to the Maxwell Garnett estimate of 
the HCM relative permittivity Chcm/o based on spherical particles of component medium a 
embedded in the host component medium b. Similarly, HS^j is equivalent to the Maxwell 
Garnett estimate of the HCM relative permittivity 6^°^^^^ based on spherical particles of 
component medium h embedded in the host component medium a. The estimate e^^i^i/^ is 
valid for /„ < 0.3, whereas the estimate e^cM/b valid for fb < 0.3; but see the footnote in 
Section 1. 

To gain insights into the asymptotic behaviour of the Wiener and Hashin-Shtrikman 
bounds, let us again restrict attention to the case of nondissipative component mediums 
(i.e., 6(15 £ 11^) • From (7), we see that Wp remains finite for all values of 5, but niay 
become infinite for 5 < since 

\Wa\^oo as S^-fy. (9) 

Jb 

In a similar vein, from (8) we find that 

\HSa\ ^ oo as 5^ ^^^^ ; (10) 

Jb 

thus, for all values of 5 < — 2 there exists a value of fb € (0, 1) at which HSa is unbounded. 
Analogously, 

\HSp\ ^ oo as 6 ^ Z^" ; (11) 

Ja "J 

so we can always find a value of /„ e (0, 1) at which HSj3 is unbounded, provided that 
<5e(-i,0). 



4 



3 Numerical results 



Let us now present, calculated values of the HCM relative permittivity e^cM) along with the 
corresponding values of the bounds Wa^p and HSa,/3, for some representative examples. Both 
nondissipative and dissipative HCMs are considered for 6 — ±3. 

3.1 Nondissipative component mediums 

The effects of dissipation may be very clearly appreciated through first considering the 
idealized situation wherein the components mediums are nondissipative [27]. Furthermore, 
although the absence of dissipation is unphysical due to the dictates of causality [28], weak 
dissipation in a particular spectral regime is definitely possible and is then often ignored [9, 
Sec.2.5]. 

Thus, it is instructive to begin with the commonplace scenario wherein both > and 
£(, > 0. For example, let €„ = 6 and ej, — 2. In Figure 1, e^^^ is plotted against /„, along with 
the corresponding Wiener bounds Wa,d and Hashin-Shtrikman bounds HSa^p- The latter 
bounds are stricter than the former bounds in the sense that 

< HS„ < e^cM < HSp < Wp. (12) 

The close agreement between ^^^^ and the lower Hashin-Shtrikman bound HS^ at low 
volume fractions /„ is indicative of the fact that HSa = Similarly, e^'^^ agrees closely 

with the upper Hashin-Shtrikman bound HSp at high values of /„ since HS^ = encM/b- 

A markedly different situation develops if the real-valued €„ and e?, have opposite signs. 
For example, the values of e^^^j^j calculated for = —6 and £6 = 2 are graphed against fa 
in Figure 2, together with the corresponding Wiener and Hashin-Shtrikman bounds. The 
Bruggeman estimate e^^M complex-valued with nonzero imaginary part for /„ < 0.82. 
This estimate is not physically reasonable. The Bruggeman homogenization formalism — 
unlike the SPFT which is its natural generalization — has no mechanism for taking coherent 
scattering losses into account. Furthermore, no account has been taken in the Bruggeman 
equation (2) for the finite size of the particles [11, 29, 30]. Therefore, the Bruggeman estimate 
of the HCM relative permittivity is required to be real-valued if the component mediums 
are nondissipative. 

While e|cM Figure 2 is complex-valued, the Wiener bounds and the Hashin-Shtrikman 
bounds are both real-valued. In accordance with (9), we see that \Wa\ — > oo as /a |. 
Similarly, \HSa\ — > oo in the limit /a ^ ;|, as may be anticipated from (10). Furthermore, 
since HSa = ^^cm/a^ Maxwell Garnett formalism is clearly inappropriate here. We also 
observe that the inequalities (12) which hold for 5 > 0, do not hold for 5 < 0. 

3.2 Weakly dissipative component mediums 

Let us now investigate e^^^ and its associated bounds when the component mediums are 
dissipative; i.e., Ca^b £ C. We begin with those cases for which S > 0: for example, we take 
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e„ = 6 + 0.3i and €6 = 2 + 0.2i. In Fi gure 3, e^Qj^ is plotted against /„, and the associated 
Wiener bounds Wa^p and the Hashin-Shtrikman bounds HSa^p are also presented. The 
behaviour of the real parts of e^^j^, Wa^js and HSa^p closely resembles that displayed in the 
nondissipative example of Figure 1. In fact, the following generalization of (12) holds: 

Re (Wa) < Re (HSa) < Re (e^'^^) < Re (HSp) < Re (Wp) . (13) 

However, this ordering (13) does not extend to the imaginary parts of e^^,^, Wa,i3 and HSa,j3- 
Turning to the cases for 5 < 0, we let €„ = — 6 + 0.3i and = 2 + 0.2i, for example. The 
corresponding Bruggeman estimate e^^^^j is graphed as function of /„, along with the Wiener 
bounds VFq-,/3 and the Hashin-Shtrikman bounds HSa^p in Figure 4. Since Im (e^ ;,) 7^ 0, the 
real parts of and HS^ remain finite, unlike in the corresponding nondissipative scenario 
presented in Figure 2. 

However, the real and imaginary parts of and HS^ exhibit strong resonances in the 
vicinity of /„ = | (for Wa) and /„ = \ (for HSa). These resonances become considerably 
more pronounced if the degree of dissipation exhibited by the component mediums is reduced. 
For example, in Figure 5 the graphs corresponding to Figure 4 are reproduced for = 
-6 + 0.003i and = 2 + 0.002i. We observe in particular that Im (e^^'^M) > 1 ioi 0.05 < 
fa ^ 0.8. Thus, the Bruggeman estimate e^^M vastly exceeds both the Wiener bounds Wa,i3 
and the Hashin-Shtrikman bounds HSa,p for a wide range of fa- Since Im {eafi) < 0.003, the 
estimates of Im(e^'^(^j) are clearly unreasonable. 

Furthermore, since the real and imaginary parts of HSa = ^ucM/a exhibit sharp resonances 
at /a = we may infer that the Maxwell Garnett formalism is inapplicable for 6 < 0. 

3.3 Highly dissipative component mediums 

On comparing Figures 4 and 5, we conclude that the Bruggeman formalism, the Weiner 
bounds and the Hashin-Shtrikman bounds become increasing inappropriate as the degree of 
dissipation decreases towards zero. This means that all three could be applicable rather well 
when the dissipation is not weak. 

Therefore, let us examine the scenario wherein the real and imaginary parts of the rel- 
ative permittivities of the component medium are of the same order of magnitude; i.e., we 
take ta = —6 + 3i and €(, = 2 + 2i. The corresponding plots of the Bruggeman estimate e^c^ 
together with the Wiener bounds Wa^p and the Hashin-Shtrikman bounds HS^^p are pre- 
sented in Figure 6. The real and imaginary parts of the Bruggeman estimate are physically 
plausible, and both lie within the Hashin-Shtrikman bounds. The Hashin-Shtrikman bounds 
themselves do not exhibit resonances, and the Weiner bounds do not exhibit strong reso- 
nances. Accordingly, we conclude that many previously published results are not erroneous, 
but caution is still advised. 
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4 Discussion 



The Bruggeman homogenization formalism is well-established in the context of isotropic 
dielectric HCMs, as well as more generally [2] . However, this formahsm was shown in Sec- 
tion 3.2 to be inapplicable for HCMs which arise from two isotropic dielectric component 
mediums, characterized by relative permittivities and e^, with 

(i) Re (e^) and Re (cb) having opposite signs; and 

(ii) |Re(€a,6) I > |Im(ea,6) |. 

Since the Bruggeman formalism provides the comparison medium which underpins the SPFT, 
it may be inferred that the SPFT is likewise not applicable to the scenarios of (i) with (ii). 

It is also demonstrated in Section 3.2 that both the Wiener bounds and the Hashin- 
Shtrikman bounds can exhibit strong resonances when the component mediums are charac- 
terized by (i) with (ii). In the vicinity of resonances, these bounds clearly do not constitute 
tight bounds on the HCM relative permittivity. As a direct consequence, the Maxwell Gar- 
nett homogenization formalism, like the Bruggeman homogenization formalism, is inappli- 
cable to the scenarios of (i) with (ii). This limitation also extends to the recently developed 
incremental [31] and differential [32] variants of the Maxwell Garnett formalism. 

If the component mediums are sufficiently dissipative then the Bruggeman formalism and 
the Hashin-Shtrikman bounds (and therefore also the Maxwell Garnett formalism) provide 
physically plausible estimates, despite the real parts of the component medium relative 
permittivities having opposite signs — as shown in Section 3.3. The explicit delineation of 
the appropriate parameter range (s) for the Bruggeman formalism and the Hashin-Shtrikman 
bounds is a matter for future investigation. 

Bounds can, of course, be violated by a formalism if the underlying conditions for the 
formalism are in conflict with those used for deriving the bounds. Sihvola [22] has catalogued 
the following conflicts: 

(a) Bounds derived for nondissipative component mediums can be invalid for the real parts 
of either or e^cM/fe ^ composite medium containing dissipative component 
mediums. 

(b) Percolation cannot be cannot be captured by the Maxwell Garnett formalism [5, 33]. 
Hence, the Hashin-Shtrikman bounds, being based on the Maxwell Garnett formalism, 
can be violated by the Bruggeman estimate e^cM fo^ ^ percolative composite medium. 

(c) The derivations of bounds generally assume that the particles in a composite medium 
have simple shapes. If the particle shapes are complicated, the composite medium may 
display properties not characteristic of the either of the component mediums. For in- 
stance, magnetic properties can be displayed when the particles in a composite medium 
have complex shapes [34, 35] , even though the component mediums are nonmagnetic. 
Clearly, the magnetic analogs of W^, Wp, HSa and HSp are then inapphcable. 
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(d) Wa, Wp, HSa, and HSp as well as their magnetic analogs are also invalid prima facie 
when the component mediums exhibit magnetoelectric properties [3, 11, 36]. 

(e) Bounds derived for electrically small particles become inapplicable with increasing 
frequency, due to the emergence of finite-size effects [29]. Even the concept of homog- 
enization becomes questionable with increasing electrical size [2, p. xiii]. 

In contrast, the bounds and the homogenization formalisms studied in this paper share 
the same premises; yet, a conflict arises in certain situations because the bounds exhibit 
resonance while the homogenization estimates do not. 

5 Concluding remarks 

As several conventional approaches to homogenization are not appropriate to the HCMs 
arising from component mediums characterized by (i) with (ii), there is a requirement for 
new theoretical techniques to treat this case. This requirement is all the more pressing, given 
the growing scientific and technological importance of new types of metamaterials [1, 26]. 

Acknowledgement. We thank two anonymous reviewers for comments that led to the 
improvement of this paper. 
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0.2 0.4 0.6 0.8 1 

fa 

Figure 1: The Bmggeman estimate e^^^j (solid line) plotted against /„ for = 6 and e?, = 2. 
Also plotted are the Wiener bounds, (thick dashed line) and Wg (thin dashed line), 
and the Hashin-Shtrikman bounds, HSa (thick broken dashed line) and HSp (thin broken 
dashed hne). 
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Figure 2: The real (top) and imaginary (bottom) parts of the Bruggeman estimate e^^^j 
(sohd hne) plotted against /„ for = —6 and = 2. Also plotted are the real parts of 
the Wiener bounds, Wa (thick dashed line) and W/j (thin dashed line), and the Hashin- 
Shtrikman bounds, HS^ (thick broken dashed line) and HSp (thin broken dashed line). The 
imaginary parts of W^^p and HSa,i3 are null-valued. 
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Figure 3: The real (top) and imaginary (bottom) parts of the Bruggeman estimate e^cM (sohd 
hne) plotted against fa for ea = 6+0. 3i and eb = 2+0. 2i. Also plotted are the Wiener bounds, 
Wa (thick dashed line) and Wfs (thin dashed line), and the Hashin-Shtrikman bounds, HSa 
(thick broken dashed line) and HSp (thin broken dashed line). 
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Figure 4: As Figure 3 but for e„ = -6 + 0.3i and e;, = 2 + 0.2i. 
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Figure 5: As Figure 3 but for e„ = -6 + 0.003i and = 2 + 0.002i. 
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